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Abstract: We perform a non-perturbative study of pure gauge theory in a two dimensional 
non-commutative (NC) space. On the lattice, it is equivalent to the twisted Eguchi-Kawai 
model, which we simulated at N ranging from 25 to 515. We observe a clear large- A scaling 
for the 1- and 2-point function of Wilson loops, as well as the 2-point function of Polyakov 
lines. The 2-point functions agree with a universal wave function renormalization. Based on 
a Morita equivalence, the large- A double scaling limit corresponds to the continuum limit 
of NC gauge theory, so the observed large-A" scaling demonstrates the non-perturbative 
renormalizability of this NC field theory. The area law for the Wilson loops holds at small 
physical area as in commutative 2d planar gauge theory, but at large areas we find an 
oscillating behavior instead. In that regime the phase of the Wilson loop grows linearly 
with the area. This agrees with the Aharonov-Bohm effect in the presence of a constant 
magnetic field, identified with the inverse non-commutativity parameter. 
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1. Introduction 

The concept of spaces with NC coordinates has been known for a number of decades (see, 
e.g. (I]]), before the mathematicians worked out a precise formulation of field theory on such 
geometries In the recent years, these NC field theories experienced a tremendous boom, 
which reflects a wide-spread belief that this variation of standard (commutative) field theory 
is well motivated. The main reason for this popularity is that they are obtained as low 
energy limits from string theories and M-theory with antisymmetric tensor fields ||, |B| . 
Moreover, there may be practical applications to the quantum Hall effect ||. For recent 
reviews we refer to refs. 0]. 

As an important deviation from standard field theory, the non-commutativity tensor 

©/ii/ = i[x v ,Xfj] (x^ : hermitean operators) 

= 0e M „ (hid = 2) (1.1) 

introduces some non-locality inside a short range, which is characterized by the scale of the 
parameter 8. This non-locality implies deep conceptual problems, but on the other hand 
it is considered a virtue which could be crucial for the link to quantum gravity as well as 
string and M-theory. Moreover, there was originally some hope that this property could 
help to weaken or even remove the UV divergences |IJ. 

However, it is exactly the question of renormalizability where NC field theory runs 
into problems which are even worse than in standard field theory. UV singularities are still 
present, and in addition new IR singularities are mixed in (this is discussed for instance in 
ref. ||). In perturbation theory, this notorious UV/IR mixing affects certain non-planar 
diagrams || . Intuitively one may grasp from eq. ( |1 . 1[ ) that short distance effects in one di- 
rection translate into long range effects in the other directions and vice versa, which causes 
this kind of mixed singularities in non-planar diagrams, which do not occur in standard 
field theory. This property poses severe difficulties for the renormalization of perturba- 
tion theory, which are not overcome yet. For instance, scalar fields can become unstable 



(tachyons) due to IR effects as the non-commutativity is switched on |J. The attempts to 
renormalize perturbation theory include the methods known from standard field theory (see 
e.g. ref. [10]), the application of Wilson's renormalization group technique fljTfl , controlling 



IR divergences in the framework of supersymmetry [12| and the Hartree method [13|. In 
spite of some plausibility arguments in favor of perturbative renormalizability, it is an open 
question if NC quantum field theories do really have finite UV and IR limits. There are 
even claims that basic NC field theories, like NC QED, are not renormalizable [14]. 



Our approach is the non-perturbative study of NC field theory by means of Monte Carlo 
simulations. In this paper we present — as a pilot project — numeric results for 2d NC 
gauge theory. 1 In this particular model, the UV/IR mixing does not occur perturbatively 
since there are no UV divergent diagrams in the commutative case. However, such mixing 
effects are still possible on the non-perturbative level, and indeed our study is going to reveal 
their existence. For interesting analytical work in this field, see refs. [|l6f . In particular, 
Paniak and Szabo worked out an explicit form for the corresponding partition function. 
Due to Morita equivalence, 2d NC U(n) gauge theory can be identified with the twisted 
Eguchi-Kawai model (TEK) |17| in the large-TV limit Q. The TEK at finite N corresponds 
to lattice NC gauge theory |ll| , as we briefly review in section |2|. In section [| we present our 
simulation results for the TEK, which do reveal that 2d NC gauge theory is renormalizable. 
We also observe the similarities and differences from commutative gauge theory and discuss 
their interpretation. Our conclusions are given in section [|. 

2. The twisted Eguchi-Kawai model as a formulation of NC gauge theory 

Back in 1982, Eguchi and Kawai conjectured that standard U(iV) and SU(iV) lattice gauge 
theories may be equivalent to their dimensional reductions to d = (one point) in the 
large- iV limit [|20| . Then all link variables are replaced by U Xill — ► U^, and the plaquette 
action simplifies to 

Sek[U) = -N/3 TrfaUyUpfy , (2.1) 

where fj,, v = 1, . . . , d. Eguchi and Kawai proved that the Schwinger-Dyson equations - 
and therefore the Wilson loops — are unaltered under dimensional reduction, assuming 
that the U(l) d symmetry of the phases (or center symmetry) is not spontaneously broken. 
This is correct in d = 2 of course, but in d > 2 it only holds at strong coupling [21|, so 
the proof for the equivalence to the reduced model (EK equivalence) has only a limited 
application. However, it was shown recently that a variant of the 4d Eguchi-Kawai model 
(EK model) surprisingly obeys the area law in a significant range of scale, so EK equivalence 
may still be in business in some finite regime [E^j. 

In order to avoid the problem with the limited viability of the proof for EK equivalence, 
the EK model was soon modified such that the phase symmetries are preserved. This was 



achieved by the "quenched EK model" |21j| , and in even dimensions also by the TEK |1 

Stek[U] = -N0 Yl Z ^ TrfaUvUlU?) . (2.2) 



1 While our simulations were running, a related project was suggested by S. Profumo TIE 



The factor Zn V represents the twist, 



where A; is an integer and L = N 2 l d is the size of the system. The usual formulation sets 
k = 1. Further variants of the TEK have been invented |23|, |24l]. For numerical studies in 



d > 2, see refs. g7|, |2§] (d = 4) and ref. |2§ (d = 6). 

Here we are concerned with the 2d TEK. It seemed for a long time that adding a 
twist is not highly motivated in d = 2, since there even the ordinary EK model coincides 
with lattice gauge theory in the planar limit, which was solved analytically by Gross and 



Witten [27]. This agreement was shown already in the famous ref. f20|| , and it seemed likely 
to hold also with twist at N — > oo. Still the 2d model was simulated with the k = 1 twist 
first in ref. [p8] and later in ref. [29], where the finite N effects were studied; but they did 



still not provide much motivation for including the twist. 2 

However the situation changed suddenly due to a new interpretation of the TEK as 
an equivalent description of NC gauge theory. The latter can be formulated if we multiply 
the fields by star products, 

S = \J d 2 xTr(F^(x)*F^(x)), 
Ffu, = d^A v - d u A^ + ig(A^ *A U -A U * . (2.4) 



For completeness we mention that the star product is defined as 



1 d d 

2 ^ v dx^ dy v 



f(x)-kg{x) = exp ( -O^-ET^— ) f( x )9(y) 



(2.5) 



x=y 



and that under a star gauge transformation 

A^x) -» V(x) * A^x) * V{x)l - -V{x) * -£-V{xf , 

9 ox^ 

V{x)*V{xy = 1, (2.6) 



the lagrangian Tr F^ v -kF^ u is not invariant, but the action in eq. fl2.4p is star gauge invariant. 

This new interpretation of twisted reduced models holds in the stringent sense of 
"Morita equivalence". Roughly speaking, this means that the algebras involved in two 
models can be mapped onto each other one-to-one, in a way which maps each module and 
its structure properly from one algebra to the other. 3 

This kind of equivalence was found between NC gauge theory and twisted large-A^ 
reduced models (T^|. That relation was established by embedding the (dynamically gener- 
ated) coordinates and momenta of the reduced model into matrices. These matrices can be 
mapped on functions, where the trace turns into an integral and the star products arise, so 



2 A very recent reconsideration of the corresponding EK equivalence is presented in ref. |30[ . 

3 As an even stronger relation, "gauge Morita equivalence" implies that also the connections between 



two modules are mapped onto each other |3l| 



that one arrives at action ( |2.4[) . The construction of NC U(n) gauge theories (for certain 
n € {1, 2, . . .}) works out in this way at N = oo. 4 Ref. |33|] pointed out that this transition 
corresponds to a Morita equivalence. At finite N the conditions cannot be matched at the 
boundaries. 



However, in refs. IS] this mapping was refined and it turned out that certain reduced 
twisted U(iV) models at finite N correspond to NC Yang-Mills theories on the lattice. The 
corresponding lattice formulation itself was also elaborated in refs. [19]. 5 From the UV/IR 



mixing it can be understood that such a formulation appears necessarily on a periodic lat- 
tice, so that the singularities are consistently regularized, and also the non-commutativity 
parameter 8 takes discrete values. The continuum limit of lattice NC gauge theory coin- 
cides precisely with the large- iV limit of the reduced (twisted, commutative) model. This 
observation completes the explicit non-perturbative ("constructive") formulation of NC 
gauge theory, in a form which can be studied by numerical simulations. This would be 
very difficult in the direct lattice formulation, due to the problems related to the star pro- 
duct: first it is complicated to introduce star unitary link variables and a measure for them, 
then one would still have to find a way to compute a discrete version of the star products 
in action (p.4j) . 

The U(iV) symmetry of the twisted reduced models corresponds to the gauge invari- 
ance ( |2.6[ ) in the NC gauge theories. In the twisted reduced model the formulation of U(iV) 
invariant observables (Wilson loops, in particular) is straightforward, and the mapping on 
NC gauge theories led also to the formulation of NC Wilson loops |52| . This formulation 
works perfectly on the lattice as well Various aspects of gauge invariant observables 



have been investigated in refs. [35|. 



In order to remain in the framework of the discrete Morita equivalence of refs. [19] 
one has to take the large-iV limit from specific finite TV" TEK formulations. For the 2d 
TEK this prescription is different from ref. |29], which used k = 1 and N even. For the 
equivalence to lattice NC gauge theory with rank-one gauge group, the twist in eq. (f2.2; ) 
has to be taken as 

Z 12 = Z* 2X = exp Li^±^\ , (2.7) 
which implies 2k = L + 1, hence L must be odd (and L = N in two dimensions). Then 



there exists an explicit map from matrices onto NC fields [36]. Corresponding to N x N 
matrices, one obtains fields on a NC torus, which is discretized into a N x N periodic lattice. 
The number of degrees of freedom matches exactly. The non-commutativity is encoded in 
the discretized version of the star product, which is used to define products of fields in 
the action. If we introduce a lattice spacing a, the (dimensionful) non-commutativity 
parameter takes the form 

8 = -Na 2 . (2.8) 

In the context of EK equivalence, one takes the limit N — > oo at fixed /?, which ensures 
the equivalence of the TEK model to U(oo) commutative lattice gauge theory. This is 

4 On the other hand, there are severe problems in formulating SU(n) gauge theory on a NC space. 
5 For earlier work on a lattice formulation of NC gauge theory, see e.g. ref. [p4|. 



called the planar large- iV limit, since only planar diagrams survive. Then the continuum 
limit should be taken by sending a — > along with (3 — > oo in a particular way dictated 
by the coupling constant renormalization. In the present d = 2 case, the coupling constant 
scaling is purely canonical (see below) 

P(xa~ 2 . (2.9) 

Now interpreting the TEK as a lattice NC gauge theory, the planar limit corresponds to 
9 = oo due to eq. ( |2.8[ ). This is a non-perturbative account of the well-known equivalence 
between 9 = oo NC field theory and the corresponding (commutative) planar field theory. 

In this paper, we are interested in a continuum limit of the lattice NC gauge theory at 
finite 9. According to eq. ( |2.8[ ), this requires the large- N limit to be taken simultaneously 
with the continuum limit a — > 0. 6 This also implies that non-planar diagrams are included, 
but in this specific model they do not involve perturbative UV/IR mixing, as we mentioned 
in the introduction. As we are going to explain in the next section, [3 has to be scaled as 
a function of the lattice spacing a in exactly the same way as in the planar theory. Hence 
in d = 2 we are going to search for a double scaling limit keeping the ratio N//3 constant. 
The question of renormalizability of 2d NC gauge theory at finite 9 can be answered by 
studying this large- N limit of the TEK. If the observables converge to finite values, we can 
conclude that 2d lattice NC gauge theory does have a finite continuum limit. Our results 
will also be presented in the next section. 



A double scaling limit of the untwisted EK model was searched for in ref. [29]. Large- N 
scaling of observables was indeed observed at a fixed ratio N//3, but this combination of 
parameters was identified empirically. The TEK, on the other hand, was studied only at 
fixed f3 in ref. and did not reveal systematic finite N effects. However, we repeat that 
the TEK was studied there with twist parameter k = 1 and even values for N. 

3. Numeric results 

In this section we present our numeric simulation results for the 2d TEK at odd N, with 



the twist specified in eq. (2/7). In the form of eq. (|2.2j ), the action cannot be simulated with 



a heat bath algorithm because it is not linear in the matrices U^. Following ref. [37] we 



introduce an auxiliary matrix field Q, which enters in a gaussian form and which linearizes 
the action in the U^, 

Stek[U, Q}=NpY^ Tr [QUQw ~ QU + tvvPJJy) - 

jl<V 

't:Mul + t*uM)Qj . (3.1) 



The auxiliary field Q consists of general complex N x N matrices, and their integration 



reproduces the TEK as given in eq. (2^2), if = Z^ v . In this respect, and in all other 



technical points of the simulation, we follow the methods which were used and described 



6 The physical extent of the NC torus Na goes to infinity in this limit. In order to obtain a finite torus 
in the continuum limit, one would have to use a different formulation of lattice NC gauge theory MM. 



N 


25 


35 


55 


85 


125 


195 


255 


515 


# of configurations 


9000 


8400 


20000 


8000 


3400 


260 


120 


18 



Table 1: Overview of the statistics in our simulation of the TEK. 



in ref. [29|. For instance, we also update SU(2) subgroups of the matrices, in the spirit 



of Cabbibo and Marinari [38]. 



Table [l| displays the number of configurations that we generated at various N. Note 
that the few configurations at our largest values of N provide conclusive signals only for 
part of the observables that we evaluated, see below. In most of our study we fixed 

N 

- = 32, (3.2) 

which corresponds to 

= -^-— = ~w2.55 (3.3) 

47T [3 IT 

in the continuum limit. 

As a point of orientation for the interpretation of our data, we consider the exact 
solution of planar U(oo) and SU(oo) lattice gauge theory in d = 2 [27]. There the Wilson 
loop of a rectangle I x J follows the area law 

w(I x J) = exp(-K(fi)IJ) , (3.4) 

where the (dimensionless) string tension k((3) is given by 

«(0H i\ (3-5) 




4/3 j y - 2- 
In particular, the internal energy amounts to 

E(P)= W (lxl) = \ P \ , (3.6) 

[ 1 4/3 P — 2 

which reveals a third order phase transition. 7 Eq. ( |3.4| ) shows how one should tune the 
bare coupling constant f3 as a function of a when one takes the continuum limit a — ► 0. In 
fact in this case the scaling is exact if one takes the lattice spacing a to be 



a = v U(P). (3.7) 

Then the Wilson loop becomes an a-independent function of the physical area a 2 IJ. From 
eq. ( |3.7| ) we see that j3 should be sent to infinity in the continuum limit as in relation ( |2.6| ) . 

The EK equivalence implies that the exact result Q3.4| ) is reproduced by considering the 
corresponding Wilson loops in the TEK in the large- iV limit at fixed (3. More specifically, 



we define 'Wilson loops' in the TEK model as [17 



W, U (I x J) = {Z liV ) I - J Tr(uj 1 U J v Ul I Ui J ) , (3.8) 



7 Note that we are always in the /3 > 1/2 regime for the values of iV that we are going to consider. 



where /, J = 1,2, For / = J = 1 it reduces to the plaquette, which appears in the 

action ( |2.2| ). Then the EK equivalence states that 

lim h\V 12 {I x J)) = u,(/ x J) , (3.9) 

JV— >oo iv 

at fixed /3. Here the term w(I x J) is the expectation value of Wilson loop in the U(oo) 
lattice gauge theory, given by eq. ( |3.4| ), 

Let us now consider Wilson loops in NC gauge theory. Since the Wilson loop ( p.g| ) 
in the TEK model is U(iV) invariant, it defines a star gauge invariant observable upon 



mapping the matrices onto NC fields [32|. This is precisely what we consider as the Wilson 
loop in NC gauge theory on the lattice jOJ. Now we take the large- N limit together with 
the continuum limit a — > in order to obtain a finite non-commutativity parameter 6 in 



eq. (^8). We define a "lattice spacing" as in the planar theory, i.e. as in eq. (3/7). Then 
the Wilson loop still agrees with the planar result (^OJ) when the physical extent of the 
Wilson loop is much smaller than the non-commutativity scale yd. (In the limit 8 — > oo it 
agrees with the planar result for arbitrarily large area, as we mentioned before.) Therefore, 
the scaling in the small area regime requires the lattice spacing to be chosen in the same 
way as in the planar theory. The first interesting question is now whether the scaling of 
the Wilson loop extends also to larger areas. 

In the EK model it was observed that square shaped Wilson loops converge faster to 
the known exact results in the planar limit than other rectangles with the same area p9| . 
Hence we also focus on square shaped Wilson loops here, which corresponds to the I = J 
case in definition ( |3.8[) . We define the normalized Wilson loop 

W(I) = ^(W 12 (IxI)), (3.10) 

which should agree with w(I x I) in eq. ( |3.4D when al <C \/9. Due to the presence of the 
twist, there is no invariance under the parity transformation U\ — ► U2, U% — ► U±. As a 
consequence, is complex in general, and Wu = W21, hence the Wilson loop depends 
on its orientation. The real part represents an average over both orientations. Figure [j] 
shows the real part Re[W(J)] as a function of the physical loop area. 

Large- N scaling is clearly confirmed. At small up to moderate areas, the Wilson loop 
follows the area law of planar lattice gauge theory given in eq. ( |3.4p , but at large areas it 
deviates and the real part oscillates around zero instead. This observation is consistent over 
a wide range of N. Remarkably, not even the absolute value decays monotonously; at large 
areas it seems to fluctuate around an approximately constant value. This is illustrated in 
figure || (on top). 

Figure || (below) shows that beyond the Gross- Witten regime, the phase increases 
linearly (in the physical area) in the large- N limit. Additional measurements at N//3 = 16, 
24 and 48, corresponding to different values of 6, reveal that the phase $ of the Wilson 
loop W(I) (beyond the Gross- Witten regime) is, in fact, given to a very high precision by 
the simple relation 

$ = 4, (3-H) 
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Figure 1: The real part of the normalized Wilson loop W(I) against the physical area for various 
values of N. At small areas Re[W(/)] follows the Gross- Witten area law, but at larger areas it 
oscillates around zero. The dashed lines are plotted to guide the eye. 



where A is the physical loop area, see figure ||. This has been confirmed also for other 
(non-square) rectangular Wilson loops, which shows that the effect does not depend on 
the shape of the Wilson loop. 8 Indeed the formula ( 3.1l| ) agrees with the Aharonov- 
Bohm effect in the presence of a constant magnetic field B = 1/8 across the plane. This 
is reminiscent of the description of NC gauge theory by Seiberg and Witten ||. Our 
numeric results seem to support a picture of that kind, which is also known in solid state 
physics f|. 

The observed large area behavior of the Wilson loop confirms that the continuum limit 
of NC gauge theory is different from any ordinary (commutative) gauge theory, hence we 
have found a new universality class. Naively one might think that in the large area regime 
A ^> y/0, the effect of non-commutativity is invisible. The fact that we do observe the 
contrary can be understood as a manifestation of non-perturbative UV/IR mixing. This 
observation also implies that the limit 9 — > does not commute with the continuum limit. 

In the double scaling limit of the untwisted EK model, the expectation value of the 
Wilson loop is real, and it remains positive even at large physical areas [29]. This means 



8 Note, however, that the expectation values of Wilson loops with the same area but with different shapes 
have in general different absolute values in the double scaling limit. 
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Figure 2: The polar coordinates of the complex Wilson loop W(I) plotted against the physical 
area A = a 2 1 2 . At small areas it is approximately real and follows the Gross- Witten area law (solid 
line), c.f. figure |l]. At larger areas the absolute value does not seem to decay any more, and the 
phase agrees with the formula ^ = A/6, c.f. figure |[ 



that the two models, twisted and untwisted, yield qualitatively different double scaling 
limits, although they become identical in the planar large- N limit. 
Figure § shows the connected Wilson loop 2-point function 



= (W 12 (I x I)W 2 i(I x /)) - (W 12 (I x I))(W 21 (I xI))£R, (3.12) 

1, . . . ,N. Here we include a wave function 

(3.13) 



again plotted against the physical area, for / 
renormalization factor 



The exponent —0.6 was found to be optimal for G 2 to scale. Indeed it leads to a neat 
large- N scaling in some regime. 

Next we consider the Polyakov line 



PJl) = Tr(u'), 



(3.14) 



which is also U(iV) invariant and therefore has an interpretation as a star gauge invariant 
observable in NC gauge theory. In fact Polyakov lines are mapped to open Wilson lines 



i 1 1 1 r 




5 10 15 20 25 



a 2 I 2 

Figure 3: The phase of the Wilson loop W(I) is plotted as a function of the physical area A — a 2 I 2 
for N = 125 and various values of N/(3, corresponding to 9 — 1.27, 1.91,2.55,3.82. For large areas 
(A > 0(9)), the phase <& agrees with the formula ^ — A/9, represented by the straight lines. 



which carry non-zero total momentum [32]. Their momentum p is related to the separation 
vector v between the two ends of the line. In general the relation is given by = Q^ u p u 
modulo the periodicity of the torus. In the present case, the Polyakov line P^I) corre- 
sponds to a momentum mode with p v = 2ir£/ (Na), where the integer £ is given by 1/2 and 
by (1 + N)/2 for even and odd /, respectively. In what follows, we plot the results against 
al for even / = 2, . . . , N — 1. 

The phase symmetry 9 makes (P^(I)} vanish, but the connected n-point functions 
(n > 1) of Polyakov lines are sensible observables. In figure || we show the 2-point function 10 

Gi P \l) = (P,(I)P^(I)) , (3.15) 

where P-^(I) = Tr^U^ 1 ). Note that there is no disconnected part in Again we insert 

the wave function renormalization which was optimal for the Wilson 2-point function, 

G% ] -» /T - 6 Gf ] . (3.16) 

As a function of the physical length al, the result is consistent with large- scaling, as well 
as a universal wave function renormalization. A similar wave function renormalization was 
also observed in the EK model 2S], where the optimal factor in relation ( 3.16| ) is modified 



to /r - 65 . 

Our statistics was not sufficient to extract conclusive signals for the connected parts 
of 3-point functions or even 4-point functions. For completeness we present nevertheless 



9 In the terminology of NC gauge theory, this corresponds to momentum conservation. 
10 The choice of the direction fj, is irrelevant. In practice we average over both possibilities in order to 
enhance the statistics. 




0.001 




Figure 4: The 2-point function G% defined in eq. (3.12), with the renormalization factor 

given in eq. (3.13) (bottom) and without (top), plotted against the physical area, for I = 1, . . . , N. 
The solid line connects the results for N = 195. 



our result for 

Gf\l) = {P,(I)P,{I)P^(2I)) (3.17) 

(there is still no disconnected part). In figure || we show our data after re-scaling, if we 
assume a universal wave function renormalization, which is not ruled out. This means that 
any n-point function is renormalized as G n — > (3~°' 6n / 2 G n (for n > 2). 
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Figure 5: The 2-point function (I), defined in eq. ( 3.15| ), with the renormalization factor given 
in eq. ( 3.13| ) (bottom) and without (top), plotted against the physical distance, for I = 2, . . . , N — 1. 
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4. Conclusions 

We have investigated the 2d TEK in the form which allows for a mapping onto 2d NC 
gauge theory on the lattice, by discrete Morita equivalence. In this representation as a 
twisted reduced model, NC gauge theory takes a form which can be simulated numerically. 
It requires the specific twist given in eq. ( [2.7| ) and odd values of N. 
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Figure 6: The 3-point function G 3 (I), defined in eq. fl3.17|) and rescaled by a universal renor- 
malization factor, plotted against the physical distance. 

The first conclusion from our simulation results is that we do observe a double scaling 
limit as N, (3 — > oo. This corresponds to the continuum limit of the NC rank-one gauge 
theory, which has therefore also finite observables. This observation demonstrates the 
non-perturbative renormalizability of 2d NC rank-one gauge theory. 

The Wilson loop follows an area law at small physical areas, and in this regime NC 
gauge theory agrees with planar standard gauge theory. However, at larger areas the Wilson 
loop becomes complex and the real part (the mean values over both loop orientations) 
begins to oscillate around zero. The phase is proportional to the physical area enclosed by 
the Wilson loop, irrespectively of its shape, and the coefficient of proportionality is given by 
the inverse of the non-commutativity parameter 9. This agrees with the Aharonov-Bohm 
effect in the presence of a constant magnetic background field 



(4.1) 



Our results support this law, which is also a key element of the Seiberg and Witten de- 
scription of NC gauge theory ||. Moreover, the same law also occurs in condensed matter 
physics: if one assumes a plane crossed by a constant magnetic flux, then the electrons in 
this plane can be projected to the lowest Landau level in a NC space, where 9 = hc/eB Q. 

The behavior of the Wilson loop at large areas implies that we have found a qual- 
itatively new universality class. At first sight, it may look surprising that the non- 
commutativity — which introduces a short-ranged non-locality in the action — changes 
the IR behavior of the gauge theory completely. However, this effect does not appear un- 
natural in the presence of UV/IR mixing. It is remarkable that this mixing effect does not 
occur in the perturbative expansion of this model, hence our results for large Wilson loops 
reveal a purely non-perturbative UV/IR mixing. 
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For the connected Wilson loop 2-point function, as well as the 2-point function of the 
Polyakov line, we can confirm the large- N scaling. These observables are in agreement with 
a universal wave function renormalization, which yields a factor /3 -0 6 for a connected 2- 
point function. This is similar to — but not identical with — the universal renormalization 
factor which was found for the (untwisted) EK model. 11 

Finally we should mention that in d > 2, NC field theories are often formulated such 
that the time coordinate is not included in the non-commutativity. Including it causes 
especially severe conceptual problems regarding causality etc. Therefore it would be im- 
portant to confirm non-perturbative renormalizability also in higher dimensions, where 
non-commutativity is introduced only between spatial coordinates. We are currently sim- 
ulating the 3d NC 4 model [40|, which can also be mapped onto a twisted reduced matrix 
model [32|. A self consistent one-loop analysis in d = 4 suggests the appearance of an 



interesting phase diagram [13]. 
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